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Abstract 

Integration of Hamiltonian systems by reduction to action-angle variables has 
proven to be a successful approach. However, when the solution depends on el¬ 
liptic functions the transformation to action-angle variables may need to remain 
in implicit form. This is exactly the case of the simple pendulum, where in order 
to make explicit the transformation to action-angle variables one needs to resort 
to nontrivial expansions of special functions and series reversion. Alternatively, 
it is shown that the explicit expansion of the transformation to action-angle vari¬ 
ables can be constructed directly, and that this direct construction leads naturally 
to the Lie transforms method, in this way avoiding the intricacies related to the 
traditional expansion of elliptic functions. 
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1. Introduction 

The simple gravity pendulum, a massless rod with a fixed end and a mass 
attached to the free end, is one of the simplest integrable models. However, it pro¬ 
vides a useful didactic system that may be available in most students laboratories 
to be used in different levels of physics [1], Besides, this physical model serves 
as the basis for investigating many different phenomena which exhibit a variety 
of motions including chaos [2], and has the basic form that arises in resonance 
problems [3]. The case of small oscillations about the stable equilibrium position 
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is customarily studied with linearized dynamics, and was for a long time the basis 
for implementing traditional timekeeping devices. 

The dynamical system is only of one degree of freedom, but the motion may 
evolve in different regimes, and one must resort to the use of special functions to 
express its general solution in closed form [4]. In fact, this dynamical model is 
commonly used to introduce the Jacobi elliptic functions [5]. 

The traditional integration provides the period of the motion as a function of 
the pendulum’s length and the initial angle [6]. Alternatively, the solution can 
be computed by Hamiltonian reduction, a case in which action-angle variables 
play a relevant role [7]. In particular, they are customarily accepted as the cor¬ 
rect variables for finding approximate solutions of almost integrable problems by 
perturbation methods [8]. 

When the closed-form solution of an integrable problem is expressed in terms 
of standard functions, the transformation to action-angle variables can be made 
explicitly in closed form, as, for instance, in the case of the harmonic oscillator. 
However, if the solution relies on special functions, whose evaluation will depend 
on one or more parameters in addition to the function’s argument, the action-angle 
variables approach may provide the closed-form solution in an implicit form. This 
fact does not cause trouble when evaluating the solution of the integrable problem, 
but may deprive this solution of some physical insight. Besides, in usual perturba¬ 
tion methods the disturbing function must be expressed in the action-angle vari¬ 
ables of the integrable problem, thus making necessary to expand the (implicit) 
transformation to action-angle variables as a Fourier series in the argument of the 
special functions. These kinds of expansions are not trivial at all, and finding them 
may be regarded as a notable achievement [9, 11, 10]. 

In the case of elliptic functions, the normal way of proceeding is to replace 
them by their definitions in terms of Jacobi theta functions, which in turn are 
replaced by their usual Fourier series expansion in trigonometric functions of the 
elliptic argument, with coefficients that are powers of the elliptic nome [12, 13]. 
This laborious procedure is greatly complicated when the modulus of the elliptic 
function remains as an implicit function of the action-angle variables, a case that 
requires an additional expansion and the series reversion of the resulting power 
series. 

On the other hand, the expansion of the transformation to action-angle vari¬ 
ables can be constructed directly. Indeed, the assumption that the transformation 
equations of the solution are given by Taylor series expansions, with the require¬ 
ment that this transformation reduces the pendulum Hamiltonian to a function of 
only the momenta, and the constrain that the transformation be canonical leads 


2 


naturally to the Lie transforms procedure [14, 15]. The procedure is illustrated 
here for the rotation regime of the simple pendulum [16, 17]. Indeed, after rear¬ 
ranging the Lie transforms solutions as a Fourier series, we check that both ki nds 
of expansions match term for term. Therefore, there is no need of making any 
numerical experiment to show the convergence of the Lie transforms solution for 
the pendulum’s rotation regime [17]. 

In the case of the pendulum’s oscillation regime the series expansion of the 
closed form solution as a Fourier series in the action-angle variables still can be 
done. However, these expansions provide a relation between an angle whose os¬ 
cillations are constrained to a maximum elongation, and an angle that rotates, for 
this reason not fulfilling the condition required by the Lie transforms method that 
original and transformed variables are the same when the small parameter van¬ 
ishes. On the contrary, the Lie transforms solution for the oscillation regime is 
commonly approached after reformulating the pendulum Hamiltonian as a per¬ 
turbed harmonic oscillator by making use of a preliminary change of variables 
to Poincare canonical variables. This case is well documented in the literature 
[19, 16] and is no tackled here. 

2. Hamiltonian reduction of the simple pendulum 

The Lagrangian of a simple pendulum of mass m and length / under the only 
action of the gravity acceleration g is written L-T-V where, noting 0 the angle 
with respect to the vertical direction, V = mgl{ 1 - cos 0) is the potential energy 
and T = (l/2)/or is the kinetic energy, where I = ml 1 is the pendulum’s moment 
of inertia and co = 0, where the over dot means derivation with respect to time. 

The conjugate momentum to 6 is given by 

0 = ^7 = ml 2 6 
80 

That is, 0 = Ico is the angular momentum. 

Hence, the usual construction of the Hamiltonian 77 = 0 0(0) - X(0,0(0)) 
gives 

^ = \mP + mgl ^ ~~ C ° S ^ 
which represents the total energy for given initial conditions 77(0 0 , 0o) = E. De¬ 
pending on the energy value the pendulum may evolve in three different regimes: 

• 0 < E < 2 mgl, the oscillation regime, with a fixed point of the elliptic type 
at E = 0 (0 = 0, 0 = 0). 
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• E = 2mgl => 0 = ±2ml 2 -\[gjl cos(0/2), the separatrix, with fixed points of 
the hyperbolic type 0 = 0, 0 = ±n. 

• E > 2 mgl, the rotation regime 


2.1. Phase space 

For each energy manifold *H = E, the Hamiltonian (1) has a geometric inter¬ 
pretation as a parabolic cylinder P E - {(x,y, z) 

x = 1 sin 6 , y = / cos 6 , 

and 

2 E 

h =- 21 = const. 

mg 

Besides, the constraint x 2 + y 1 = / 2 makes that the phase space of the sim¬ 
ple pendulum is realized by the intersection of parabolic cylinders, given by the 
different energy levels of the Hamiltonian (1), with the surface of the cylinder 
C = {( x,y,z ) | v 2 + y 2 = l 2 } of radius /. This geometric interpretation of the phase 
space is illustrated in Fig. 1. 

Alternatively, typical trajectories on the cylinder are displayed by means of 
simple contour plots of the Hamiltonian (1), as illustrated in Fig. 2 where the 
trajectories are traveled from left to right for positive heights (0 > 0) and from 
right to left for negative heights (0 < 0). 

The traditional solution to the motion of the pendulum is approached by the 
direct integration of Hamilton equations 


| Z 2 - 2y = h}, with 
0 


m/p 1 / 2 ' 


. m © . m 

0 = -W = - rafesin( '’ 

which are usually written as a single, second order differential equation 


0 + y sin 0 = 0. (2) 

Note that the summand mgl in Eq. (1) is a constant term that does not affect 
the dynamics derived from Hamilton equations. Hence, this term is sometimes 
neglected, which implies a trivial displacement of the energy by a constant level. 
Besides, the Hamiltonian can be scaled by the pendulum’s moment of inertia to 
show that it only depends on a relevant parameter. Nevertheless, the dimensional 
parameters are maintained in following derivations, because, in addition to the 
physical insight that they provide, the simple and immediate test of checking di¬ 
mensions is very useful in verifying the correctness of the mathematical develop¬ 
ments. 
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Figure 1: Geometric construction of the pendulum’s phase space. Parabolic cylinders Ve = 
{(x,y,z)\z 2 -2 y - const.}. 
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Figure 2: Phase space of the simple pendulum (q = yflm 2 gP). 
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2.2. Hamiltonian reduction 

Alternatively to the classical integration of Eq. (2), the flow can be integrated 
by Hamiltonian reduction, finding a transformation of variables 

T :(0',0')->(0,0), (3) 

such that the Hamiltonian (1) when expressed in the new variables is only a func¬ 
tion of the new momentum. Namely, 

T \ r H = f H(0(6', 0'), 0(0', 0')) = ®(- 0'). (4) 


The solution of Hamilton equations in the new variables 


O' 


5® 

S0 7 ’ 


fv 50 
0 =- 

06 ' 


0, 


is trivial 

0' = const, 6' = 6' 0 + co't, (5) 

where the frequency a>' = <9®/<90' only depends on 0' and, therefore, is constant. 
Plugging Eq. (5) into the transformation T in Eq. (3) will give the time solution 
in the original phase space (0,0). 

The Hamiltonian reduction in Eq. (4) can be achieved by the Hamilton-Jacobi 
method [ 1 8] , in which the canonical transformation T is derived from a generating 
function S = S{6, 0') in mixed variables, the “old” coordinate and the “new” 
momentum, such that the transformation is given by 


e = — . 
30' 00 


( 6 ) 


The functional expression of 0, given by the last of Eq. (6) is replaced in the 
Hamiltonian (1) to give the partial differential equation 

i(§) +2mgl sin2(9/2) = o<0 ' ) ’ <7) 


where the trigonometric identity 1 - cos 6 = 2 sin 2 (0/2) has been used. 

The Hamilton-Jacobi equation (7) is a partial differential equation that can be 
solved by quadrature 


S = V2 ml 2 J' ^®(0') - 2 mgl sin 2 (0/2) d 0. 


( 8 ) 
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But, in fact, the generating function S does not need to be computed because the 
transformation (6) only requires the partial derivatives of Eq. (8). Thus, 


* = r 

d® Jo 


d§ 


2^0(0') -2mg/sin 2 (tf/2) 


0 = 


0(0') - 2 mgl sin (0/2), 


(9) 

( 10 ) 


where d = 0 corresponds to the upper extreme of the square root in Eq. (9). 

Remarkably, the quadrature in Eq. (9) can be solved without need of choosing 
the form of the new Hamiltonian in advance [17, 20]. 

2.3. Rotation regime: 0(0') > 2 mgl 
Equation (9) is reorganized as 

where the partial differentiation has been replaced by the total derivative in view 
of the single dependency of O on 0', and 


F(0 


rt> 

,k 2 ) = 

Jo 


d6 


Vl - k 2 sin 2 6 

is the incomplete elliptic integral of the first kind of amplitude 

and elliptic modulus 


formation in mixed variables 

O' = 


<p = 0/2, 

(id 

J 2mgl <i. 

(12) 

y 0(0') 

of the simple pendulum is 

solved by the trans- 

d<D 

lc — F(0,*r), 
d0' 

(13) 

: VfW‘-* 2sin >' 

(14) 
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which is, in fact, a whole family of transformations parameterized by ®. For 
instance, a “simplifying” option could be choosing ® = 0 ,2 /(8 ml 2 ), which results 
in 



0 


2 yjm 2 gP 


sn (8',k 2 ), 

j dn(0', k 2 ), 
k 


(15) 

(16) 


where k = 4 ^Jm 2 gP /&, from Eq. (12). Note that f = am (O', k 2 ), where the Jacobi 
amplitude am is the inverse function of the elliptic integral of the first kind, and sn 
and dn are the Jacobi sine amplitude and delta amplitude functions, respectively. 

From the definition of k in Eq. (12), one may express the reduced Flamiltonian 
in the standard form 

® = 2 mgl^. (17) 

Then, by simple differentiation, 


d® 

d0' 


„ ,1 d k 

-ding l— 


k 3 d0'' 


(18) 


which turns Eq. (13) into 


6' = -4 mgl 



d k 
g k 2 d0 


;F(M 2 ). 


(19) 


3. Action-angle variables 

Note that, in general, the new variables ( 0',®') defined by the mixed trans¬ 
formation in Eqs. (13) and (14) may be of a different nature than the original 
variables (0,0): an angle and angular momentum, respectively. Indeed, the new 
Hamiltonian choice leading to Eqs. ( 1 5)— ( 1 6) assigns 6' the dimension of time. 
However, choosing a transformation that preserves the dimensions of the new 
variables (0', 0') as angle and angular momentum —the so-called transformation 
to action-angle variables— may provide a deeper geometrical insight of the so¬ 
lution and, furthermore, is specially useful in perturbation theory. Alternative 
derivations to the one given below can be found in the literature [3]. 

Hence, the requirement that 0' = 0'(0,0) be an angle, given by the condition 
§d0' = In [8], is imposed to Eq. (13). That is, when 0 completes a period along 
an energy manifold 'Hid, 0) = E then O' must vary between 0' () and 0' 0 + In. 
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Note in Eq. (11) that when 9 evolves between 0 and 2/r, (p = <p(9 ) evolves 
between 0 and n. Therefore, using Eq. (13), the angle condition reads 



where 


F(0, k 2 ) = 0, F(tt, k 2 ) = 2F(tt/2, k 2 ) = 2K (k 2 ). 


Hence, 



( 20 ) 


which is replaced in Eq. (13) to give 



( 21 ) 


On the other hand, by eliminating dO/d0' between Eqs. (18) and (20) results 
in separation of variables. Then, 0' is solved by quadrature to give 



( 22 ) 


where E(k 2 ) = E(7t/2, k 2 ) is the complete elliptic integral of the second kind, and 



is the incomplete elliptic integral of the second ki nd. 

It deserves noting that the elliptic modulus cannot be solved explicitly from 
Eq. (22), and hence the reduced Hamiltonian O must remain as an implicit func¬ 
tion of 0' in the standard form of Eq. (17). However, the rotation frequency of 6' 
is trivially derived from Hamilton equations by using Eq. (20), namely 



d ff dO 


n 


(24) 



l kK(k 2 )' 


In summary, starting from ( 9 ,0), corresponding action-angle variables (O', 0') 
are computed from the algorithm: 

1. compute E = El id, 0) from Eq. (1); 
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2. compute k from Eq. (12) with O = E; 

3. compute 0' from Eq. (22), and 0' from Eq. (21), where 0 = 0/ 2. 

On the other hand, starting from ( 0 ',&), the inverse transformation is com¬ 
puted as follows. 

1. solve the implicit equation (22) for k\ 

2. invert the elliptic integral of the first kind in Eq. (21) to compute 

0 = am(7r _1 K(fc 2 ) 0', k 2 }; 

3. 6 = 20, and 0 is evaluated from Eq. (14). 


4. Series expansions 


In spite of the evaluation of elliptic functions is standard these days, using 
approximate expressions based on trigonometric functions will provide a higher 
insight into the nature of the solution and may be accurate enough for different 
applications. 

Thus, Eqs. (21) and (14) are written 


0 


0 


2 am(w, k 2 ), 


2ml 2 



dn (u,k 2 ). 


(25) 

(26) 


which is formally 
for the argument 


the same transformation as the one in Eqs. (15) and (16), except 


u = F(0, k 2 ) = 



n 

which now fulfils the requirement that 0' be an angle. 

The elliptic functions in Eqs. (25) and (26) are expanded using standard rela¬ 
tions, 3 


dn (u, k 2 ) 
am (u, k 2 ) 


n 


1 


+ 


2 Zt 

n= 1 


q 


+ q 


m 2 ) 

+ 2 y q 

2 ^ n( 1 + q 2n ) 

n-\ ± 


2n 


cos (n0') 


sin (n0'). 


(27) 

(28) 


3 See http://dlmf.nist.gov/22.llE3 and http://dlmf.nist.gov/22.16.E9 for 
Eqs. (27) and (28), respectively. 
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where the elliptic nome q is defined as 

q = exp [-7 tK( 1 - k 2 )/K(k 2 )], (29) 

and k must be written explicitly in terms of 0'. This is done by expanding 0' as 
given by Eq. (22) in power series of k, viz. 


0' 


2 yjgPm 
k 

441 

"65536 




4851 

1048576 


k n + 0(k 13 ) 


J1 5 _ k * 

16384 


Then, by series reversion, 


k = 2 


1 


f+ 4 6 


-e 3 + 


161 

~64' 


239 

~64 


e 5 + <9(e 6 ) 


(30) 


where the notation 


m 




e- • <31) 

has been introduced. Needles to say that e must be small in order to the expansions 
make sense. 

Once the required operations have been carried out, the transformation to 
action-angle variables given by Eqs. (25) and (26) is written as the expansion 


0 


, . 11 , 247 

^ + |1 + T6 e+ 256 eUsmfl 


+ -^-e 2 |e 3 sin 36' 
48 64 


+ (i + ii f2 ) f2sin2< '' + 

+ _l_ e 4 S i n4 S - + _l_ f5s i n5e ' + o ( A 


(32) 



r i 2 

15 4 

( t 3 o 

39 4 \ 

0 = 0' 

1 - -e 2 - 

-e 4 + 

1 + — e" + 

6 4 

2 

32 

l 16 

256 / 


1 1 


+ I t + t e 2 6 2 cos 20' + I — + — e 1 1 e 3 cos 38' 


4 4 


1 


2 1 ,3 


16 64 


+ J-e 4 cos 40' + cos 50' + 0(e 6 ) 
64 256 


(33) 


Proceeding analogously, the standard Hamiltonian in Eq. (17) is expanded like 


<D = 


0 


/2 


2 l 2 m 


! 1 2 5 

1 + 2 f + 32 


f 4 + 


64 




(34) 
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5. Action-angle variables expansions by Lie transforms 


In spite of modem algebraic manipulators provide definite help in handling 
series expansion and reversion, the procedure described in the previous section 
can be quite awkward. Then, it emerges the question if the expanded solution 
in Eqs. (32), (33) and (34) could be constructed directly. The answer is in the 
affirmative, and dealing with elliptic functions and series reversion can be totally 
avoided by setting the rotation regime of the pendulum as a perturbed rotor. In¬ 
deed, since previous expansions are based in the convergence of Eq. (30), the 
function e = 6(0') must be small enough. Hence, neglecting the constant sum¬ 
mand mgl, Eq. (1) can be written as 


H = 


©2 

2ml 2 



2 m 2 gl 3 

0 2 


cos 6 


(35) 


and in those cases in which 2 m 2 gP/& 2 <sc 1, that is 


mgl 


0 2 

2 ml 2 


the reduction of the Hamiltonian (35) can be achieved as follows. 


5.1. Deprit’s triangle 

First, Eq. (35) is written 


'H = 


0 2 

2 ml 2 


a: cos 6, 


where the physical parameter k = mgl is small. 4 Then, the pendulum Hamiltonian 
can be written as the Taylor series 


with the coefficients 


'H = V —H ni) (Q,<d), 

(36) 

©2 

H °’° 2ml 2 ’ 

(37) 

H t o = - cos 6, 

(38) 

H IU o = 0 {n > 1). 

(39) 


4 In fact, one can always choose the units so that © 2 /(2ml 2 ) is of order 1 and k <s 1. 
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The reasons for using a double subindex notation will be apparent soon. 
The desired transformation T is formally written 

6 = 0(0',0';k), 

0 = 0(0\ ©';*), 


which is assumed to be analytic. Then, when this transformation is applied to 
Eq. (35), one gets 

T> = ‘H(0(0',0';k),0(0',0 , ;k);/O, 
which is expanded as a Taylor series in the new variables 


n> 0 


, f^0,n ? 
n\ 


H, 


d' ,v P 


0 ,n 


d K n 


K= 0 


(40) 


where the required derivatives with respect to the small parameter k are computed 
by the chain rule 

d^ _ m de m d0 &h 

ck “ ~50~ck + 50 ck + ~ 8 k 

We want the transformation to be canonical and explicit, and hence impose 
that T is derived from some generating function TV = TV(0,0; k). In particular, 
it is imposed that T be a Lie transform, so that it is defined by the solution of the 
differential system [19, 21] 


d6 _ <9TV d0 _ <9TV 

ck _ "50"’ ~dK~~~df' 

for the initial conditions k = 0, 


(42) 


do = 6(6', Q'- 0 ) = e\ 

0o = 0(0', 0';O) = 0'. 


The existence of the solution of Eq. (42) in a neighborhood of the origin, is guar¬ 
anteed by the theorem of existence and uniqueness of the ordinary differential 
equations. Besides, since the transformation is computed as the solution of a 
Hamiltonian system, its canonicity is also guaranteed. 

Then, taking into account Eq. (42), Eq. (41) is rewritten as 


d'F 5-K5TV 57T5TV d<H &H 

d K de 50 50 de dK x ’ 1 dK ’ 


(43) 
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where {*#, TV] is the Poisson bracket of 'H and TV. 

Replacing TV in Eq. (43) by its Taylor series expansion 


TV 


K n 

Erf* 


n> 0 


n +1 


( 0 , 0 ), 


( 44 ) 


where the subindex n +1 is now used for convenience in subsequents derivations of 
the procedure, after straightforward manipulations, one arrives to the recurrence 

Hn, q+ 1 = H n+hq + 2 r) {#„_,-*; W i+1 }. (45) 

0 <i<n ' ' 


Finally, since the terms H( )n needed in Eq. (40) are evaluated at k = 0, the 
original variables are replaced by corresponding prime variables. Full details on 
the derivation of Eq. (45), may be found in Deprit’s original reference [15] or in 
modern textbooks in celestial mechanics. 

The recurrence in Eq. (45) is customarily known as Deprit’s triangle. Indeed, 
the computation of the term H 0 n requires the computation of all intermediate 
terms in a “triangle” of “vertices” // 0 ,o, #o,«, and H n{) . For instance, for n = 4 
the recurrence in Eq. (45) results in a triangle made of: 

#o,o 

#i,o #o,i 

# 2,0 H 11 

H 3,0 #2,1 #1,2 #0,3 


# 4,0 # 3,1 # 2,2 # 1,3 # 0,4 

Deprit’s triangle is not constrained to the case of one degree of freedom Hamil¬ 
tonians and generally applies for the transformation of any function F = Fix, X; k) 
where x e R", X e R" are coordinates and their conjugate momenta, respectively. 
In particular, when the generation function is given by its Taylor series expansion, 
it is easily checked that the solution of the differential system (42) is equivalent to 
applying Deprit’s triangle to the functions 


o = V — 0„,o(0,0), 0 = V. — 0„,o(0,0) 

n\ n\ 

ri> 0 n >0 


( 46 ) 


where 0o,o = 0, 0o,o = 0, and d nfi = 0„, o = 0 for n > 0. 
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5.2. Perturbations by Lie transforms 

Deprit’s triangle provides an efficient and systematic way of computing higher 
orders of a transformation when the generating function is known. However, it is 
precisely the generating function that is needed to be computed with the only 
condition that the transformation reduces the pendulum Hamiltonian to a func¬ 
tion of only the (new) momentum. Therefore, a new Hamiltonian fulfilling this 
requirement must be also constructed. The whole procedure of constructing the 
action-angle Hamiltonian and computing the corresponding generating function 
can be done in an iterative way as follows. 

The construction of ® is an iterative procedure based on Deprit’s triangle. At 
each order n Eq. (45) is rearranged in the form 

{W n ,H 0 ,o} + H iln =H n , (47) 

which is dubbed as the homological equation, where 

• H n is known from previous computations 

• h 0m is chosen in agreement with some simplification criterion such that 

• W n is obtained as a solution of the partial differential equation (47) 

Commonly, the simplification criterion for choosing H( )n is eliminating cyclic 
variables in order to reduce the number of degrees of freedom of the new, trans¬ 
formed Hamiltonian after truncation to 0(K n ). This is exactly the case required for 
the reduction of the pendulum Hamiltonian. 

In practice, the procedure starts by reformulating all the functions involved in 
the homological equations in the prime variables, so that the new Hamiltonian is 
directly obtained in the correct set of variables. 


5.3. The pendulum as a perturbed rotor 

We start from Eqs. (36)-(39). For any term n of the generating function, 


(S)' BW 

The first order of the homological equation is 


0' BW { 
ml 2 86' 


= H ( u-H u 


(48) 
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with H\ = H l 0 . Then, H ()] is chosen by averaging terms depending on the cyclic 
variable O' 

1 r 2n ~ 

// 0 ,i = ;- H,dff = 0, 

2tt Jo 

and Wi is solved from Eq. (48) by quadrature 

ml 2 C ~ , ml 2 , 

w= or J(H 0 ,i-H 1 )dtr = -— W i<r. 

Once W\ is obtained, the first order of the transformation is calculated by 
application of Deprit’s triangle to Eq. (46) with x = 6 and X = 0. It is obtained 

0o,i = = 

K 

0 O1 = {©', VTJ = 0'-cos0', 

K 

where the abbreviation e = nrgl 3 /®' 2 is used, in agreement with Eq. (31). 

The second order of the homological equation is 


0' dW 2 

ml 2 d6 


= Hog ~ H 2 , 


(49) 


where H 2 is computed from Deprit’s triangle 

H 2 = {H xfi ,Wi) = ^sin 2 0', 
0'- 

H 0 2 is chosen by averaging 


H, 


0,2 


i r 2n 

2tt Jo 


Ho d0' = 


ml 2 

20^ 


and W 2 is solved from Eq. (49) by quadrature 

72 


Wo 


ml 1 r 

= 0 7 J 


(Ho, 2 -H 2 )d6' = 


0' 2 6 2 
2ml 2 k 2 ’ 


m 2 / 4 

40' 3 


sin 26'. 


Then, a new application of Deprit’s triangle to the variables provides the second 
order of the transformation 

1 6 2 

00,2 = {0o.uWi) + {0,W 2 ) = -—Sin 20', 

4 /c- 

00,2 = {0o,i,Wi} + {0 , ,^2} = 0 , ^^cos20 , -lj, 
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New iterations of the procedure lead to the third order 


# 0,3 

# 0,3 


0, 

e 3 / 33 . . 1 . „ . 

— — sin 6 + - sin 3# 

k 3 \ 8 8 


0i 


0,3 


= 0 ' 


e 3 19 


k 3 \8 


cos O' + - cos 36 '), 


the fourth order, 


#< 

#< 

0i 


0' 2 15 e 4 


0,4 - 

2ml 2 


e 4 19 

* 0,4 = 

K 4 \2 



* 0,4 = 

0'-| 
K 4 ' 


sin 26' + — sin 46' I, 

0' ^7 (6 cos 26' + ^ cos 46' - ^ |, 


the fifth order, 

H, 


0,5 

# 0,5 


= 0, 

e 5 / 3705 


45 3 

. . sin 6' + — sin 36' + — sin 56' I, 
k 5 \ 32 8 32 1 


0, 


0,5 


= 0 ' 


e 5 /585 


105 


15 


. . „ cos 6' + —— cos 36' + — cos 5#' 
k 5 32 8 32 


and so on. 

Now, it is simple to check that the new Hamiltonian % = 2ji>o(k"/ w 0#o,/ 
matches Eq. (34). Besides, rearranging the perturbation series as Fourier series 

# = J]^0 0 ,n(0',&) = J]soA®')sin(nd'), 

n >0 ’ n >0 

0 = ,„(#',0') = 2c o ,„(0') cos («#'), 

rc>0 * n >0 

it is simple to check that they match Eqs. (32) and (33), respectively. 


6. Conclusion 

Expanding the transformation to action-angle variables of the simple pendu¬ 
lum solution as a Fourier series, is a cumbersome task because of the elliptic 
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integrals and functions on which the solution relies upon, on the one hand, and 
the implicit character of the transformation, on the other. Alternatively, these 
Fourier series expansions can be computed by the Lie transforms method. The 
latter is straightforward, avoids dealing with special functions, provides the trans¬ 
formation series explicitly without need of series reversion, and is amenable of 
automatic programming by machine, thus allowing one to compute higher orders 
of the solution. 
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